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Abstract 

The  asymptotic  stress  and  velocity  fields  of  a  crack  propagating  steadily  and  quasi- 
statically  into  an  elastic-plastic  material  are  presented.  The  material  is  characterised  by  J2 

flow  theory  with  linear  strain-hardening.  The  possibility  of  reloading  on  the  crack  flanks  is 
taken  into  account.  The  cases  of  anti-plane  strain  (mode  III),  plane  strain  (modes  I  and  II), 
and  plane  stress  (modes  I  and  II)  are  considered.  Numerical  results  are  given  for  the 
strength  of  the  singularity  and  for  the  distribution  of  the  stress  and  velocity  fields  in  the 
plastic  loading,  elastic  unloading  and  plastic  reloading  regions,  as  functions  of  the  strain¬ 
hardening  parameter.  An  attempt  is  made  to  make  a  connection  with  the  perfectly-plastic 
solutions  in  the  limit  of  vanishing  strain-hardening. 

1.  Introduction 

Knowledge  of  the  stress  and  deformation  fields  near  the  tip  of  an  advancing  crack  is 
central  to  the  continued  development  of  engineering  fracture  mechanics,  and  important 
advances  have  been  made  in  recent  years.  Under  the  assumption  of  quasi-static  growth  and 
small-strain  ideal  elasto-plasticity,  rigorous  asymptotic  solutions  have  been  found  in  anti¬ 
plane  strain  by  Chitaley  and  McClintock  (1971);  in  plane  strain  by  Slepyan  (1974),  Gao 
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(1980),  Rice  etal.  (1980),  and  Drugan  etal .  (1982);  and  in  plane-stress  mode  II  by  Ponte 
(1985a).  Tn  addition,  a  general  discussion  of  the  subject  is  found  in  Rice  (1982). 

The  corresponding  work  for  hardening  elastic-plastic  solids  is  at  a  less  developed 
stage.  For  linear-hardening  materials  Amazigo  and  Hutchinson  (1977)  (referred  to  as  AH  in 
the  sequel)  produced  steady-state  asymptotic  solutions  using  Mises  flow  theory  for  plane- 
strain  and  plane-stress  mode  I,  and  for  anti-plane  strain.  These  solutions,  however, 
neglected  the  possibility  of  reverse  plastic  flow  on  the  crack  flanks  which  is  known  to 
occur  in  the  elastic -perfectly  plastic  solutions  mentioned  above.  In  a  recent  publication 
Zhang  etal.  (1983)  included  reloading  (as  well  as  a  Bauschinger  effect)  in  their  solution  of 
the  plane-strain  mode  I  problem. 

Further  work  on  the  hardening  elastic-plastic  steadily  moving  crack  include 
publications  by  Slepyan  (1973)  who  studied  a  mode  III  crack  propagating  in  a  linear- 
hardening  J2-deformation  theory  material,  Lo  and  Pierce  (1981)  who  considered  a  mode  in 

crack  moving  in  a  linear-hardening  comer-theory  material,  Gao  etal.  (1983)  who 
investigated  a  mode  III  crack  propagating  in  a  power-hardening  J2-flow  theory  material, 
and  Ponte  (1985b)  who  studied  the  mode  I  crack  propagating  in  a  linear-hardening  J2- 
deformation  theory  material. 

In  this  work  we  will  include  the  possibility  of  reloading  on  the  crack  flanks  of  the 
linear-hardening  J2-flow  theory  problem.  We  will  consider  plane  strain  and  plane  stress 

(modes  I  and  II ),  and  anti-plane  strain.  Thus  we  will  verify  and  expand  the  results  of 
Zhang  et  al.  for  plane-strain  mode  I,  we  will  correct  the  results  of  AH  for  plane-stress 
mode  I  and  mode  III,  and  we  will  produce  new  results  for  plane-strain  and  plane-stress 
mode  II. 

We  propose  to  use  a  simple  scheme  that  can  be  implemented  numerically  in  a 
straightforward  manner.  Thus  we  formulate  the  problem  in  terms  of  a  system  of  first  order 
ODE's  in  the  angular  variations  of  the  non-zero  components  of  the  stress  tensor  and  the 
velocity  vector  instead  of  a  single  higher  order  equation  in  the  stress  function  as  was  done 
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by  AH.  This  procedure  has  the  advantage  that  it  can  directly  make  use  of  well-developed 
Runge-Kutta  integrators,  and  also  that  the  boundary  and  continuity  conditions  are  simplest 
in  form. 

In  our  numerical  solutions  we  will  go  to  very  low  values  of  the  hardening  parameter 
in  order  to  try  to  make  a  connection  with  the  perfectly-plastic  problem.  However,  the  work 
of  Dunayevsky  and  Achenbach  (1982)  suggests  that  the  region  of  validity  of  the  power 
singularity  solution  to  the  linear-hardening  problem  vanishes  in  the  limit  of  zero  strain¬ 
hardening.  Thus  it  is  not  clear  whether  the  limiting  solutions  to  the  linear-hardening 
problem  can  in  fact  be  related  to  the  the  perfectly-plastic  solutions.  Here,  we  will  address 
these  and  related  questions. 


2.  Anti-plane  strain  (mode  III) 


2.1  Governing  equations. 

With  reference  to  Figure  2.1,  let  Xj  (i  =  1,2,3)  be  a  Cartesian  coordinate  system  of 
fixed  orientation  travelling  with  the  crack  tip  such  that  the  x3-axis  coincides  with  the 
straight  crack  front.  Also  let  ej  be  the  unit  vector  corresponding  to  the  x;  direction. 
Similarly,  let  r,  0  be  polar  coordinates  corresponding  to  xa  (a  =  1,2)  and  er,  e0  be  the 
conesponding  unit  vectors.  The  crack  tip  moves  with  velocity  V  =  Vej  with  respect  to  the 
stationary  coordinate  system  X;.  In  our  steady-state  analysis  the  crack  tip  speed  V  is 

constant  so  that  the  material  derivative  is  given  by 

O'— V()fl  (2.1) 

The  non-zero  stresses  are  xa  =  ca3,  and  the  only  non-zero  velocity  component  is 
v3.  Because  it  turns  out  to  be  simpler  from  an  algebraic  point  of  view,  we  choose  to  use  the 
polar  components  of  the  stress  vector  x  so  that  xr  xe,  and  v3  are  the  three  variables  of  this 
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problem.  Then  the  equilibrium  equation  is 
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(rxr),r  +  xe,6  =  0  (2.2) 

For  our  bilinear  flow  theory  material  the  stress-strain  relation  is 

2GD  =  T  +  Or1  -  1)  (te7  te)  t  (2.3) 

where  D  =  (l/2)Vv3  is  the  strain-rate  vector,  T  =  %'  is  the  stress-rate  vector,  Te  =  |t|  is 
the  effective  stress,  and  fl  is  either  a  for  loading  (te‘  S  0),  or  unity  for  unloading  (te'<  0), 
where  a  =  Gt/G,  the  ratio  of  the  tangent  modulus  to  the  elastic  modulus  in  shear  (Fig.  2.1), 
is  a  measure  of  the  hardening. 

It  can  be  seen  that  (2.2)  and  (2.3)  form  a  system  of  three  first  order  PDE’s  in  the 
two  stress  components  and  the  single  velocity  component.  Even  though  some  of  the 
equations  are  nonlinear  (namely  the  constitutive  equations),  all  three  equations  are 
homogeneous  in  the  independent  variables,  which  suggests  that  we  look  for  solutions  of 
the  form 


Tr(r,0)  =  T0  y  t(e)  Is  t9(r,0)  =  T0  y2(0)  rs 

(2.4) 

v3(r,0)  =  V  (tq/G)  y3(0)  is/s 

where  £  is  nondimensionalised  with  respect  to  some  measure  of  the  plastic  zone  size,  R. 
Putting  these  into  (2.2)  and  (2.3)  we  are  left  with  a  system  of  three  first  order  ODE's  in  the 


vector  y  =  (y1,y2,  y3) 


y’(0)  =  f(0,  y;  s,  &)  (2.5) 

where  the  the  components  of  f  are  given  in  Appendix  A  and  depend  on  whether  loading  or 


unloading  occurs  through  the  parameter 
&  =  a  if<j>’  £0 

=  1  if  <}>*  <  0 


(2.6) 


where 

4>'(0)  =  r  Te'/  xe  =  -s  cos0  +  sin0  [ye'(9)  /  ye(9)l  (2.7) 

xe  =  xo  y«(®)  Is  (2.8) 

and  s  is  some  as  yet  unknown  function  of  a. 


********* 
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2.2  Boundary  conditions 

Anti-plane  symmetry  ahead  of  the  crack  tip  reduces  to  the  two  independent 
conditions 

yi(0)=y3(0)  =  0  (2.9) 

and  vanishing  traction  on  the  crack  face  requires  that 

y200  =  0  (2.10) 

Thus  we  have  three  homogeneous  boundary  conditions  to  be  satisfied  by  (2.5).  It 
follows  that  the  strength  of  the  singularity,  s,  will  be  determined  as  an  eigenvalue  of  the 
problem,  and  an  amplitude  factor  corresponding  to  the  stress  intensity  factor  in  the  elastic 
problem  will  be  left  undetermined  in  this  asymptotic  solution.  We  can  then  initially  set 

y2(0)  =  l  (2.11) 

even  though  once  we  obtain  a  solution  we  will  renormalise  it  such  that 

ye(®t)  =  1  (2.12) 

where  9,  is  the  unloading  angle  (to  be  defined  in  the  next  section). 

2.3  Loading  history  of  a  near-tip  particle. 

Consider  the  state  of  stress  of  a  particle  P  at  a  given  distance  x2  above  the  crack  line 

as  it  moves  past  the  crack  tip  (or,  equivalently,  as  the  crack  tip  moves  past  it)  as  depicted  in 
Figure  2.2a.  Note  that  we  can  describe  the  location  of  this  particle  by  giving  its  X! 
coordinate,  or,  alternatively,  by  giving  its  angular  coordinate  0  =  tan-](  x2/  Xi) . 

It  is  expected  that  the  particle  will  experience  plastic  loading  ahead  of  the  crack,  but 
then  unload  elastically  for  some  critical  value  of  0,  denoted  0t,  defining  the  location  of  the 

straight  unloading  boundary.  This  is  seen  to  be  consistent  with  the  fact  that  in  the  elastic 
solution  te  has  a  maximum  for  a  value  of  0  between  zero  and  n. 

For  a  hardening  material  (0  <  a  <  1),  a  particle  in  the  elastic  unloading  zone  retains 
the  plastic  strain  state  it  had  when  it  unloaded.  Thus  it  is  expected  that  when  the  particle  is 
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deep  into  the  wake  of  the  crack  (  0  — »  n),  its  associated  stress  state  will  approach  some 
non-vanishing  wake  level,  Tw.  Plastic  reloading  takes  place  on  the  crack  flank  if  the 
effective  stress  level,  tcw,  associated  with  the  wake  stress  state  is  greater  than  the  effective 
stress  level,  t,.1  =  Toy,^)  (  x2/  sin©!)5,  that  the  particle  had  at  unloading.  Hence,  it  is 
expected  that  for  some  values  of  a,  plastic  reloading  will  start  at  some  critical  value  of  0, 
denoted  02,  defining  the  reloading  angle. 

2.4  Continuity  conditions 

We  need  to  connect  the  solutions  in  the  different  regions  through  appropriate 
continuity  conditions.  Recalling  that  in  this  problem  we  have  a  third  order  system  in  the 
two  stress  components  and  the  single  velocity  component,  we  require  three  independent 
continuity  conditions  in  these  three  variables  across  each  unloading  and  reloading 
boundary. 

Let  [  ]  denote  the  jump  in  a  quantity  as  0  increases  infinitesimally.  Then  the  traction 


component  of  the  stress  must  be  continuous 

[xe]=0  (2.13) 

and  the  displacement  field  must  be  continuous 

[u3]  =  0  (2.14) 

Also  for  a  work-hardening  material  it  is  consistent  with  the  plastic  rule  to  assume  that  the 
plastic  strain  vector  is  continuous 

[yP1]  =  0  (2.15) 

It  follows  from  (2.14)  that 

[Yrl  =  [u3,r]  =  0  (2.16) 

This  condition,  together  with  (2.15)  yields 

[Yrel]  =  0  (2.17) 

which  in  turn  implies  that 
[Tr]  =  0 


(2.18) 


and  also  that 


[tj  =  0  (2.19) 

Similarly,  it  follows  from  (2.13)  that 

[YeelJ  »  0  (2.20) 

and  that 

[Ye]  =  0  (2.21) 

Finally,  note  that  the  velocity  is  given  by 

v3  =  -Vu3)1  =  -V(cos0  u3r-  sin0  u3e)  (2.22) 

so  that  it  follows  from  strain  continuity  that  the  velocity  field  must  also  be  continuous 

[v3]  =  0  (2.23) 

Conditions  (2.13),  (2.18),  and  (2.23)  are  the  three  independent  conditions  needed, 
and  they  take  the  convenient  form  in  terms  of  angular  variations 

[yi]  =  [y2]  «  [y3]  =  o  (2.24) 


In  order  to  close  the  mathematical  problem  we  need  to  specify  the  conditions  for 

I 

unloading  and  reloading  which  determine  the  angles  0!  and  02. 

Unloading  occurs  at  0j  when  the  effective  stress-rate  of  the  particle  vanishes,  or 
(see  (2.7))  when 

<>  (0t)  =  0  (2.25)  I 

Reloading  occurs  at  02  if  the  effective  stress  of  the  particle  regains  its  unloading 
value,  or  (see  (2.8))  if 

ye(0i)  I  (sin0i)s  -  ye(02)  /  (sin02)‘  =  0  (2.26) 

i 

2.5  Numerical  integration  \ 

We  need  to  solve  the  third  order  system  (2.5)  subject  to  the  three  boundary  ! 

conditions  (2.9)  and  (2.10),  and  the  continuity  conditions  (2.24).  We  selected  a  Runge-  i 

i 

Kutta-Vemer  fifth  and  sixth  order  scheme  to  perform  the  numerical  integration. 

\ 

\ 

f 

i 

f 

i 
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We  provide  the  initial  values  of  the  dependent  variables  y^O)  through  (2.9)  and 
(2.1 1),  and  we  need  to  guess  the  value  of  s  for  a  given  a.  We  then  integrate  equations 
(2.5)  to  find  the  values  of  y^G)  for  0  <  0  <  n,  keeping  track  of  whether  unloading  or 
reloading  occurs  in  order  to  make  use  of  the  appropriate  value  of  fl.  Once  we  have  the 
values  of  y^Tt)  we  check  to  see  whether  y2(rt)  is  zero,  and  we  iterate  in  our  guess  for  s 
until  convergence  is  achieved. 

Note,  however,  that  we  cannot  use  (2.5)  to  find  numerically  the  values  of  y^O), 

but  we  can  easily  take  the  appropriate  limits  to  find  that 
yj’(O)  =  (l-s2/a)  /  (1-s/a) 

y2'(0)  =  0  (2.27) 

y3'(0)  =  -s2/a 

Similarly,  equations  (2.5)  are  numerically  ill-conditioned  at  0  =  n,  but  we  can 
integrate  them  out  to  0  =  7t  -  e,  and  make  £  as  small  as  needed  to  obtain  accurate  results. 
Since  y2(0)  is  well  behaved  near  0  =  7t,  this  approximation  presents  no  difficulties. 

2.6  Results. 

The  results  are  summarised  in  Figures  2.3  through  2.6,  and  in  Table  2.1.  Figure 
2.3  is  a  plot  of  the  strength  of  the  singularity,  s,  versus  the  square  root  of  the  hardening 
parameter,  a.  It  is  observed,  as  pointed  out  by  AH,  that  s  — a1/2  for  small  a. 

Table  2.1  presents  values  of  s,  0j  and  02  for  values  of  a  such  that  10-6  £  a  <  1 , 
and  it  also  includes  for  comparison  the  values  of  s  and  0t  obtained  by  AH.  It  is  found  that 
reloading  occurs  for  all  values  of  a  such  that  0  <  a  <  a*  where  0.05  <  a*<  0.1  (a  more 
accurate  determination  of  a*  is  an  expensive  task  not  worth  the  investment  at  this  time). 
See  Figure  2.2b  for  a  schematic  depiction  of  the  loading  history  of  a  fixed  particle.  For 
a*<  a  <  1  our  results  agree  with  those  of  AH  as  expected.  For  0  <  a  <  a*  the  appearance 
of  the  thin  reloading  sector  affects  rather  weakly  the  values  of  s  and  0(.  Finally,  observe 
that  as  a  0  the  values  of  0t  and  02,  though  still  far  from  the  corresponding  perfectly- 

i 

i 
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plastic  values  (19.69°  and  179.63°)  of  Chitaley  and  McClintock  (1971),  do  appear  to  have 
an  appropriate  trend. 

Figure  2.4  -.hows  plots  of  the  particle  effective  stress  (fixed  x2)  for  0  <  0  <  n  for 
three  values  of  a  (0.75,  0.1,  0.001)  corresponding  to  large,  moderate  and  small  strain¬ 
hardening.  The  result  for  large  strain-hardening  is  actually  quite  close  to  the  elastic  result, 
and  no  reloading  is  necessary.  As  a  becomes  smaller,  the  stress  level  ahead  of  the  crack 
(0  <  8  <  9[)  approaches  a  flat  distribution,  and  a  tiny  reloading  zone  appears  on  the  crack 

flank. 

Figure  2.5  presents  plots  of  the  angular  (fixed  r)  variations  of  the  stresses  and 
velocity  for  large,  moderate  and  small  strain  hardening.  The  curves  corresponding  to  a  = 
0.75  are  not  very  different  from  the  elastic  curves  with  the  exception  of  xr  (and  te)  which 
near  0  =  re  suddenly  become  unbounded.  For  a  =  0.1  the  angular  variations  of  the  stresses 
and  velocity  undergo  significant  changes,  and  for  a  =  0.001  they  start  taking  the  general 
shape  of  the  perfectly-plastic  results  which  are  depicted  in  Figure  2.6.  This  is  a  remarkable 
result:  even  though  the  radial  dependence  of  the  velocity  of  the  small-a  problem  (rs/s)  does 
not  approach  that  of  the  perfectly-plastic  problem  (lnr),  the  factored-out  angular  variation 
of  the  velocity  of  the  small-a  problem  does  seem  to  approach  the  perfectly-plastic  result. 

3.  Plane  strain  (mode  I) 


3.1  Formulation. 

Here  the  non-zero  stress  components  are  cap  and  a33,  and  the  non-zero  velocity 
components  are  va.  Again  for  algebraic  convenience  we  choose  to  use  the  cylindrical 
components  of  the  stress  tensor  o  and  of  the  velocity  vector  v,  so  that  cr9,  cn,  aee,  a33, 
vr,  and  ve  are  the  six  variables  of  this  problem.  Then  equilibrium  requires 
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(roie),r+  °ee,e  +  °rs  =  0  (3-2) 

For  a  bilinear  (a  =  Et/E)  flow-theory  material  the  stress-strain  relation  is 

ED  =  (l+v)X  -  v  Tr(£)  I  +  (3/2)  (fl-J-l)  (ce 7  oe)  S  (3.3) 

where  v  is  Poisson’s  ratio,  I  =  e^is  the  identity  tensor,  D  =  (l/2)[Vv+(Vv)1]  is  the 
strain-rate  tensor,  £  =  o'  is  the  stress-rate  tensor,  S  =  o  -  (1/3)  Tr(a)  I  is  the  stress- 
deviator  tensor,  ae=  [(3/2)S:S] 1/2  is  the  effective  stress,  and  &  is  either  a  for  loading 
(ac‘  >  0),  or  unity  for  unloading  (oe‘  <  0). 

Note  that  there  are  four  non-trivial  equations  in  (3.3)  (including  the  plane  strain 
condition  D33  =  0)  which  means  that  (3.1),  (3.2),  and  (3.3)  form  a  system  of  six  first  order 

PDE's  in  the  six  dependent  variables  of  the  problem  identified  above.  As  we  found  in  mode 
HI,  all  six  equations  are  homogeneous  in  the  independent  variables  which  suggests  that  we 
look  for  solutions  of  the  form 

vr(r,0)  =  V  (t0/E)  yx(0)  is/s  ve(r,0)  =  V  (x0/E)  y2(0)  is/s 
CrfjM)  =  'Co  y3(0)  rs  ^(r,©)  *  t0  y4(0)  Is  (3.4) 

oee(r,0)  =  ^  y5(0)  Is  c?33(r,0)  =  t0  y6(0)  rs 

Putting  these  in  equations  (3.1)  through  (3.3)  we  are  left  with  a  system  of  six  first  order 
ODE’s  in  the  vector  y  =  (  yu  y2, y6) 

y  (0)  =  f(  0,  y;  s,  a,  v)  (3.5) 

where  the  components  of  f  are  given  in  Appendix  A,  and  depend  on  whether  loading  or 

unloading  occurs  through  the  parameter 

&  =  a  if  <}>'  £0 


=  1  if  <J>'  <  0 

where 

$'  =  r  ce'/ce  =  -s  cos0  +  sin0  (ye'/ye)  (3.7) 

°c  =  Toye(0)is  (3-8) 

and  s  is  some  as  yet  unknown  function  of  a  and  v. 

Mode  I  symmetry  ahead  of  the  crack  tip  requires  the  four  independent  conditions 


*  -  •  ..  *  *•  «"  >  . '»  .*•  .*«i  '•  >*•  '  -  ~  _■»  * 
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y2(0)  =  y3(0)  =  y4'(0)  =  y6'(0)  =  0  (3.9) 

and  vanishing  tractions  on  the  crack  faces  requires  that 

y3(7t)  =  y5(7t:>=  0  (3.10) 

Thus  we  have  six  homogeneous  boundary  conditions  to  be  satisfied  by  equations 
(3.5).  As  before  s  will  be  determined  as  an  eigenvalue  of  the  problem  for  given  a  and  v, 
and  we  can  arbitrarily  pick  the  amplitude  factor  of  the  angular  variations. 

The  same  arguments  used  in  the  mode  III  problem  can  be  repeated  here,  and  they 
lead  us  to  expect  a  similar  near-tip  loading  history  for  this  case. 

We  next  use  the  same  development  as  before  to  determine  the  continuity  conditions 


across  an  elastic-plastic  boundary  (unloading  or  reloading). 

The  traction  components  of  the  stress  tensor  must  be  continuous 

[a*,]  =  [cee]  =  0  (3.11) 

The  displacement  field  must  also  be  continuous 

[ur]  =  [u0]  =  0  (3.12) 

and  it  is  consistent  with  the  plastic  rule  to  assume  that  the  plastic  strain  tensor  is  continuous 
[EP‘]  =  0  (3.13) 

It  follows  from  (3.12)  and  the  plane  strain  condition  that 

[Err]  =  [e33l  =  0  (3.14) 

which  together  with  (3.11)  and  (3.13)  imply  that 

[^  =  [0331=0  (3.15) 

and  also  that 

[oj  =  0  (3.16) 

We  can  next  show  that  all  the  strain  components  are  continuous,  and  therefore  that 

I  vr]  =  [ve]  =  0  (3.17) 

Conditions  (3.11),  (3.15),  and  (3.17)  are  the  six  conditions  needed,  and  in  terms  of  the 
angular  variations  they  take  the  convenient  form 

[yil  -  ty2]  =  -  =  tyel  =  0  (3.18) 


To  these  conditions  we  must  add  an  unloading  condition  which  determines  0j 

4>-(0i)  =  O  (3.19) 

and  a  reloading  condition  which  determines  02 

ye(0t)  /  (sin©!)*  -  yc(02)  /  (sin02)5  =  0  (3.20) 

The  numerical  integration  of  equations  (3.5)  was  performed  using  the  same  scheme 
used  in  the  mode  III  problem.  In  order  to  get  started  we  need  to  provide  the  values  of  y(0), 
y’(0)  and  s  for  given  a  and  v. 

We  have  already  seen  in  (3.9)  that  y2(0)  =  y3(0)  =  0.  If  we  now  pick  the 
normalisation  y4(0)  =  1,  and  guess  the  values  of  y5(0)  and  s,  we  find  using  (3.5)  that 
y i(0)  =  -s  (y4(0)/a  -  X[y5(0)  +  y6(0)]} 

y6(0)  =  a  X[y4(0)  +  y5(0)]  (3.21) 

with  X  =  v  + 

Similarly,  since  y4'(0)  =  y6'(0)  =  0,  we  can  use  (3.5)  to  find  that 
yi’(O)  =  y5'(0)  =  0 

y2'(0)  =  -yi(0)  -  s2{y5(0)/a  -  X[y4(0)+y6(0)]}  (3.22) 

y3’(0)  =  -(l+s)y4(0)  +  y5(0) 

We  then  integrate  equations  (3.5)  to  find  the  values  of  y(0)  for  0  <  0  <  jt  keeping  track  of 
whether  unloading  or  reloading  occurs.  Given  the  values  of  y3(jt)  and  y5(7t)  we  check  to 
see  whether  they  vanish,  and  iterate  in  our  guesses  for  y5(0)  and  s  until  convergence  is 
achieved.  Once  this  is  accomplished  we  renormalise  our  solutions  such  that  ye(0j)  =  1. 
Note  that  y3  and  y5  are  well  behaved  near  9  =  tc,  so  that  even  though  (3.5)  is  numerically 
ill-conditioned  at  0  =  n,  a  procedure  analogous  to  that  used  in  the  mode  IU  case  is 
appropriate. 

3.2  Results. 

In  this  case  we  produced  results  for  v  =  1/3  for  several  values  of  a.  It  was  noted  by 
AH  that  the  dependence  of  s  and  of  the  angular  variations  of  the  stresses  on  v  is  weak. 


Figure  3.1  depicts  a  plot  of  s  versus  a1/2;  for  large  strain-hardening  it  looks  very  much  like 
the  corresponding  plot  for  anti-plane  strain,  but  as  a  gets  smaller  the  curve  reverses 
curvature,  and  s  does  not  seem  to  be  approaching  zero. 

Table  3.1  shows  values  of  s,  0!  and  02  for  0.001  <  a  <  1,  and  it  also  includes  for 
comparison  the  values  of  s  and  Qx  obtained  by  AH  ignoring  reloading.  We  first  note  that 
our  results  agree  very  well  with  those  reported  by  Zhang  el  al.  (1983).  It  is  found  that 
reloading  occurs  for  0  <  a  <  a*  where  0.2  <  a*  <  0.3  .  See  Fig.  2.2b  for  a  pictorial 
representation  of  the  loading  history.  In  contrast  with  the  mode  III  results,  reloading  here 
starts  at  higher  strain-hardening  (a*  -  0.3),  the  width  of  the  reloading  sector  is  much  larger 
(31.6°  at  a  =  0.001),  and  the  values  of  s  and  0j  are  significantly  affected  by  reloading  (s  is 
weaker  by  10%,  and  0!  is  smaller  by  14%  at  a  =  0.01).  Note  further  that  for  a  =  0.001 

the  width  of  the  elastic  sector  is  only  1.35“,  and  still  decreasing  strongly.  Finally,  we 
should  remark  that  the  values  of  0j  and  02  do  not  appear  to  be  approaching  the 
corresponding  perfectly-plastic  values  (123.1*  and  160.4*  for  v  =  0.3)  of  the  solution  of 
Drugan  et  al.  (1982). 

Figure  3.2  shows  plots  of  the  particle  (fixed  x2)  effective  stress  for  0  <  0  <  Jt  for 
large,  moderate  and  small  strain-hardening  corresponding  to  values  of  a  of  0.75,  0.1  and 
0.001,  respectively.  For  large  strain-hardening  small  changes  from  the  elastic  solutions 
occur,  and  no  reloading  is  necessary.  For  moderate  strain-hardening  well-defined  changes 
take  place:  a  smooth  unloading  (continuous  ae’),  and  a  sudden  (discontinuous  ae  ) 
reloading  are  clearly  observed.  As  we  make  a  smaller,  the  elastic  unloading  sector  virtually 
disappears,  and  even  though  the  stress  level  ahead  of  the  crack  increases,  it  is  still  far  from 
approaching  a  flat  distribution  corresponding  to  the  perfectly-plastic  yield  condition. 

Figure  3.3  presents  plots  of  the  angular  (fixed  r)  variations  of  the  stresses 
Ort,  an,  Oqq,  O33,  and  ae  for  large,  moderate  and  small  strain-hardening.  The  curves  for 

a  =  0.75  are  not  very  different  from  the  elastic  curves  except  that  aw  a33  and  ae  become 
unbounded  near  0  =  71.  Figure  3.4  shows  the  angular  variations  of  the  Cartesian 
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components  of  the  stress  for  a  =  0.001.  It  can  be  seen  that  the  angular  variations  of  the 
stresses  seem  to  be  approaching  the  general  shape  of  the  stresses  in  a  Prandtl  field.  But  for 
a  =  0.001  the  agreement  is  still  not  very  good;  for  instance  oe(0)  changes  by  50%  between 
zero  and  0!  instead  of  remaining  constant  throughout. 

Figure  3.5  shows  plots  of  the  angular  variations  of  the  Cartesian  components  of  the 
velocity  for  large,  moderate  and  small  strain-hardening  (a  =  0.75,  0.1,  0.01).  It  is 
observed  that,  contrary  to  what  happens  in  mode  Eli,  the  angular  variations  of  the  velocity 
do  not  remain  of  the  same  order  of  magnitude  as  those  of  the  stresses,  and  do  not  seem  to 
be  approaching  a  limit  as  a  -»  0  .  This  is  an  understandable  consequence  of  the  fact  that  s 
does  not  vanish  in  the  limit  in  this  case. 


4,  Plane  strain  (mode  II) 


4.1  Formulation. 

The  formulation  of  this  problem  is  identical  to  that  of  the  mode  I  problem  except  for 
the  boundary  conditions  ahead  of  the  crack  tip.  Mode  II  symmetry  requires  the  four 
independent  conditions 

yi(0)  =  y4(0)  =  y5(0)  =  y6(0)  =  0  (4.1) 

It  is  convenient  here  to  pick  the  normalisation  y3(0)  =  1,  and  to  guess  the  values  of  y2(0) 

and  s.  We  can  then  use  equations  (3.5)  to  find  that 

y,'(0)  =  (1-s)  y2(0)  -  2s2  [(1+v)  +  (3/2)(a-i-l)]  y3(0) 

y2’(0)  =  y3’(0)  =  0  (4.2) 

y5'(0)  =  — (s-t-2)  y3(0) 

plus  two  other  more  complicated  equations  for  y4'(0)  and  y6'(0)  which  we  will  not  write 
out  in  detail. 
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Given  the  values  of  y(0)  and  y'(0)  we  can  start  integrating  equations  (3.5)  making 
use  of  the  continuity  conditions  developed  in  Section  3  to  connect  regions  of  different 
physical  behavior.  As  before  we  need  to  check  whether  y3(rc)  and  y5(7t)  are  zero,  and  to 
iterate  in  the  guesses  of  s  and  y2(0)  until  convergence  is  achieved.  When  we  are  finished 
we  renormalise  the  solution  in  the  usual  way  such  that  ye(9i)  =  1. 

4.2  Results. 

In  this  case  we  produced  results  for  v  =  1/3  for  several  values  of  a.  Figure  4.1 
shows  a  plot  of  s  versus  a1/2;  it  looks  very  much  like  the  corresponding  mode  III  plot,  and 
unlike  the  mode  I  plot  it  is  seen  that  s  — a1/2  for  small  a. 

This  case  has  the  new  feature  that,  since  ye(6)  has  two  peaks  (recall  the  elastic 
result),  two  additional  regions  appear  and  the  two  angles  03  and  04  are  required  to  describe 

the  beginnings  of  the  new  reloading  zone  and  the  new  unloading  zone,  respectively.  The 
two  angles  B}  and  02  retain  their  old  meanings.  Table  4.1  shows  values  of  s,  0  lt  02,  03, 
and  04  for  0.001  <  a  <  1  .  It  is  found,  as  depicted  in  Fig.  2.2c,  that  for  a+  <  a  <  1  (0.3 
<  a+  <  0.5)  there  is  a  four  region  solution:  a  plastic  loading  zone  ahead  of  the  crack, 
followed  by  a  wide  elastic  unloading  sector,  followed  by  a  thin  plastic  reloading  sector, 
followed  by  another  elastic  unloading  sector.  For  a*<  a  <  a+  (0.05  <  a*  <  0.1)  there  is 
a  two  region  solution,  and  for  0  <  a  <  a*  reloading  occurs  on  the  crack  flank,  and  there  is 
a  three  region  solution.  It  should  be  noted  that  for  a+  <  a  <  1  two  unloading  conditions 
are  needed,  and  the  appropriate  obvious  notational  changes  were  made  in  equations  (3.19) 
and  (3.20);  for  0  <  a  <  a+  the  old  formulation  holds.  By  comparison  to  the  mode  III 
results  we  find  that  reloading  on  the  crack  flank  occurs  roughly  at  the  same  level  of 
hardening  (moderate),  the  singularity  is  a  little  stronger  and  the  width  of  the  reloading 
sector  on  the  crack  face  is  smaller.  By  comparison  to  the  mode  I  results  we  find  that  for 
large  strain-hardening  the  singularity  is  a  little  weaker,  and  reloading  starts  later  and  is  less 
important.  Finally,  we  observe  that  the  trend  of  the  values  of  0j  and  02  as  a  ->  0  are 
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compatible  with  the  corresponding  mode  II  perfectly-plastic  values  (9.57°  and  179.82')  of 
Slepyan  (1974)  for  a  material  with  a  Tresca  yield  surface  with  v  =  1/2. 

Figure  4.2  shows  plots  of  the  particle  effective  stress  as  a  function  of  0  for  large, 
moderate,  and  small  strain-hardening.  It  is  noted  that  for  large  strain-hardening  small 
changes  from  the  elastic  result  take  place,  and  the  isolated  plastic  reloading  zone  is  clearly 
observed.  As  a  becomes  smaller  this  initial  reloading  zone  disappears,  and  for  still  lower 
values  of  a  a  thin  plastic  reloading  zone  appears  on  the  crack  flank.  At  the  same  time  a  flat 
variation  in  oe  is  approached  ahead  of  the  crack. 

Figure  4.3  presents  plots  of  the  angular  variations  of  the  stresses  a^,  om  o00, 
c33,  and  ae  for  large,  moderate,  and  small  strain-hardening.  For  large  a,  the  usual 
observations  are  made.  For  small  a,  we  observe  from  this  and  Figure  4.4  (showing  the 
angular  variations  of  the  Cartesian  components  of  the  stress)  that  the  stress  distribution 
agrees  in  character  with  the  perfectly-plastic  results  of  Slepyan  (1974);  we  observe,  for 
instance,  a  centered  fan  distribution  developing  ahead  of  the  crack. 

Figure  4.5  shows  plots  of  the  angular  variations  of  the  Cartesian  components  of  the 
velocities  vt  and  v2.  We  remark  that  for  small  a  they  agree  in  form  with  the  perfectly- 

plastic  results  of  Slepyan  (1974).  They  are  of  the  same  order  of  magnitude  as  the  angular 
variations  of  the  stresses,  and  they  approach  a  constant  level  in  the  elastic  unloading,  and 
plastic  reloading  zones  as  a  ->  0. 


5.  Plane  stress  (mode  I) 


5.1  Formulation. 

The  variables  for  this  geometry  are  the  three  non-zero  components  of  the  stress 
tensor  0,9,  a„,  Oqq,  and  the  two  in-plane  components  of  the  velocity  vector  vr,  v0  (v3  does 

not  enter  the  formulation  of  the  problem).  The  governing  equations  are  the  two  equilibrium 
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equations  (3.1),  (3.2),  and  the  three  in-plane  components  of  the  stress-strain  relation  (3.3) 
where  we  set  c33  equal  to  zero.  These  equations  form  a  system  of  five  first  order  PDE's  in 

the  five  dependent  variables  identified  above.  If  we  now  put 

vr(r,0)  =  V  (T0/E)  yi(0)  ve(r,0)  =  V  (t0/E)  y2(0)  i* 

(5.1) 

<Wr«e)  =  t0  y3(0)  Is  aa(r,0)  =  x0  y4(0)  i*  ce0(r,0)  =  T0  y5(0)  Is 

into  the  five  equations  mentioned  above,  we  obtain  a  system  of  five  first  order  ODE's  in  the 
vector  y  =  (yj,  y2,  ....  ys) 

y'(0)  =  f(0,  y;  s,  a,  v)  (5.2) 

where  the  components  of  f  are  given  in  Appendix  A,  and  depend  on  whether  loading  or 
unloading  occurs  in  the  usual  way. 

Mode  I  symmetry  ahead  of  the  crack  requires  the  three  independent  conditions 

y2(0)  =  y3(0)  =  y4‘(0)  =  0  (5.3) 

and  vanishing  tractions  on  the  crack  faces  requires  that 

y3(*)  =  YsOO  -  0  <5-4) 

Thus  we  have  five  homogeneous  boundary  conditions  to  be  satisfied  by  our  fifth  order 
system.  As  usual,  s  will  be  determined  as  the  eigenvalue  of  the  problem,  and  it  will  depend 
only  on  a .  Similarly,  we  will  be  able  to  arbitrarily  pick  the  amplitude  factor  of  y(0). 

Next  we  can  use  the  usual  arguments  to  show  that  the  five  dependent  variables  of 
the  problem  are  continuous  across  any  unloading  or  reloading  boundary  so  that 

[yil  =  to!  =  -  =  [ysl  =  0  (5-5) 

Similarly,  we  use  the  unloading  condition  (3.19)  and  the  reloading  condition  (3.20)  to 
determine  0j  and  02  respectively. 

Once  again  we  use  the  same  numerical  scheme  to  perform  the  integration,  but  we 
need  to  get  started  by  providing  the  values  of  y(0),  y'(0),  and  s  for  a  given  a. 

From  (5.3)  we  have  that  y2(0)  =  y3(0)  =  0.  If  we  now  pick  y4(0)  =  1,  and  guess 
the  values  of  y5(0)  and  s,  we  find  using  (5.2)  that 
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yi(0)  =  -s  [y4(0)/a  -  X  y5(0)] 
and  given  that  y4‘(0)  =  0  we  find  that 
yi’(O)  =  y5'(0)  =  0 
y2'(0)  =  — y  i (0)  -  s2  (y5(0)/a  -  X  y4(0)]  (5.7) 

y3'(0)  =  -(l+s)  y4(0)  +  y5(0) 

We  then  integrate  equations  (5.1)  numerically  to  find  the  values  of  y(0)  for  0  <  0  <  7t. 
Given  the  values  of  y3(jt)  and  y5(7t),  we  check  to  see  whether  they  vanish,  and  if  they  do 
not,  we  iterate  in  our  guesses  for  y  5(0)  and  s  until  convergence  is  achieved.  The  solutions 
are  finally  renormalised  in  the  usual  way  (ye(0j)  =  1) . 

5.2  Results 

In  this  case  we  produced  results  for  a  value  of  v  equal  to  1/2,  and  several  values  of 
a.  AH  showed  through  their  formulation  of  the  problem  that  both  s  and  the  angular 
variations  of  the  stresses  are  independent  of  v.  Only  the  angular  variations  of  the  velocities 

depend  on  v.  Figure  5.1  shows  a  plot  of  s  versus  a1/2,  and  we  can  see  that  s - a1/2  for 

small  a. 

Table  5.1  shows  values  of  s,  0j,  and  02for  0.0001  <  a  <  1,  and  it  also  includes  for 
comparison  the  values  of  s  and  0!  obtained  by  AH  neglecting  reloading.  The  loading 
history  is  depicted  in  Figure  2.2b.  It  is  found  that  reloading  occurs  for  all  values  of  a  such 
that  0  <  a  <  a*  (0.005  <  a*<  0.01),  but  that  the  width  of  the  reloading  sector  is  so  small 
(less  than  0.001*)  that  the  values  of  s  and  0j  are  virtually  unaffected  by  reloading.  As 
expected,  for  a*  <  a  <  1  our  results  agree  very  well  with  those  of  AH. 

Figure  5.2  depicts  plots  of  the  particle  effective  stress  for  0  <  0  <  n  for  the  three 
usual  values  of  a.  It  is  noted  that  for  large  strain-hardening  a  small  perturbation  from  the 
elastic  result  occurs.  As  a  becomes  smaller  the  tiny  reloading  zone  on  the  crack  face 
appears,  and  the  stress  level  ahead  of  the  crack  approaches  a  flat  variation  corresponding  to 
the  perfectly-plastic  yield  condition. 
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Figure  5.3  presents  plots  of  the  angular  variations  of  the  stresses  for  large, 
moderate,  and  small  strain-hardening.  The  curves  for  large  strain-hardening  show  little 
variation  from  the  elastic  result  with  the  exception  of  Gn  (and  oe)  which  becomes 
unbounded  as  0  — >  K.  The  results  for  small  strain-hardening  show  that  the  stresses  ahead 
of  the  crack  are  approaching  a  centered  fan  distribution.  Figure  5.4  presents  the 
corresponding  angular  variations  of  the  Cartesian  components  of  the  stress  tensor. 

Figure  5.5  shows  the  angular  variations  of  the  Cartesian  components  of  the  velocity 
for  the  three  values  of  a.  We  observe  that  they  stay  of  the  same  order  of  magnitude  as  the 
stresses,  and  we  can  also  see  that  they  are  nearly  constant  in  the  elastic  unloading  sector. 


6.  Plane  stress  (mode  II) 


6.1  Formulation. 

The  formulation  of  this  problem  is  identical  to  that  of  the  corresponding  mode  I 
problem  with  the  exception  of  the  boundary  conditions  ahead  of  the  crack  which  for  mode 
II  become 

y,(0)  =  y4(0)  =  y5(0)  =  0  (6.1) 

It  is  convenient  here  to  pick  the  normalisation  y3(0)  =  1,  and  we  must  guess  the  values  of 
y2(0)  and  s  for  a  given  a.  We  can  then  use  equations  (5.2)  to  find 

yjXO)  =  (1-s)  yi(0)  -  2s2[(l+v)  +  (3/2)(a-i-l))  y3(0) 

y2'(0)  =  y3'(0)  =  0  (6.2) 

ys'(O)  =  -(s+2)  y3(0) 

plus  a  more  complicated  expression  for  y4'(0)-  With  these  results  we  can  start  integrating 
equations  (5.2)  numerically,  iterating  in  the  values  of  y2(0)  and  s  until  a  solution  is  found 
satisfying  y3(7t)  =  y5(7t)  =  0.  Finally,  we  renormalise  the  solution  in  the  usual  way. 
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6.2  Results. 

In  this  case  we  produced  results  for  v  =  1/2,  and  several  values  of  a.  Figure  6.1  is 
a  plot  of  s  versus  a1/2,  and  it  shows  the  usual  behavior  of  s  — a1/2  for  small  a. 

Table  6.1  shows  values  of  s,  0j,  02,  ©3,  and  04  for  0.001  <  a  <  1.  The  four 

angles  have  the  same  interpretation  as  in  the  plane-strain  mode  II  problem.  The  loading 
history  is  shown  in  Figure  2.2d.  A  four  region  solution  is  found  for  a*  <  a  <  1  (0.1  <  a* 
<  0.2),  a  five  region  solution  for  oc+  <  a  <  a*  (0.01  <  a+  <  0.05),  and  a  three  region 
solution  for  0  <  a  <  a+.  We  remark  that  for  a*  <  a  <  1  two  unloading  and  one  reloading 
conditions  are  needed,  and  for  a+  <  a  <  a*  two  unloading  and  two  reloading  conditions 
are  needed.  Thus  for  a+  <  a  <  1  the  appropriate  obvious  notation  changes  were  made  in 
equations  (3.19)  and  (3.20),  and  for  0  <  a  <  a+the  old  formulation  holds.  Reloading 
starts  at  intermediate  strain-hardening,  and  the  size  of  the  reloading  zone  is  small  but  not  as 
small  as  in  mode  I.  The  strength  of  the  singularity  is  between  that  of  the  mode  I  (plane 
stress)  and  mode  III  problems.  As  a  -» 0  the  values  of  0t  and  02  seem  to  be  approaching 

limiting  values  that  compare  favorably  with  the  perfectly-plastic  values  (13.31*  and 
179.61°)  of  Ponte  (1985a). 

Figure  6.2  depicts  plots  of  the  particle  effective  stress  as  a  function  of  0  for  large, 
moderate,  and  small  strain-hardening.  It  is  noted  that  for  large  strain-hardening  small 
changes  from  the  elastic  result  occur,  and  the  isolated  reloading  zone  is  observed.  As  a 
becomes  smaller  the  initial  reloading  zone  disappears,  and  a  new  reloading  zone  appears  on 
the  crack  flank.  Simultaneously,  the  stress  approaches  a  constant  level  in  the  plastic  sector 
ahead  of  the  crack. 

Figure  6.3  presents  plots  of  the  angular  variations  of  the  stresses  for  the  three  usual 
values  of  a.  Figure  6.4  presents  a  plot  of  the  corresponding  variations  for  the  perfectly- 
plastic  material.  It  is  seen  that  the  agreement  between  the  small-a  and  perfectly-plastic 
results  is  very  close. 
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Figure  6.5  shows  plots  of  the  angular  variations  of  the  Cartesian  components  of  the 
velocity  vector  for  the  linear-hardening  material.  Figure  6.6  shows  the  corresponding 
variations  for  the  perfectly-plastic  material.  Once  again,  remarkably,  it  appears  that  the 
linear-hardening  variations  approach  the  corresponding  perfectly-plastic  variations  as  a 
approaches  zero. 

7.  Discussion  of  the  small-a  limit 

It  is  natural  to  separate  the  perfectly-plastic  problems  into  two  groups  according  to 
their  solutions.  For  the  anti-symmetric  modes  (II  and  HI),  asymptotic  solutions  exist  with 
continuous  stress  and  velocity  fields,  and  they  have  roughly  the  same  configuration:  a 
(singular)  centered  fan  plastic  sector  ahead  of  the  crack,  followed  by  a  (non-singular) 
elastic  sector,  followed,  in  turn,  by  a  (non-singular)  constant  stress  plastic  sector  on  the 
crack  flank.  For  the  symmetric  mode  (I),  fully-continuous  solutions  with  the  above,  or 
other  configurations  have  not  been  found,  and  presumably  do  not  exist.  For  plane-strain 
mode  I,  a  solution  has  been  found  (Rice  et  al.,  1980)  which  admits  a  discontinuity  in  the 
tangential  component  of  the  velocity  vector  at  the  plastic-plastic  boundary  (0  =  45°) 
separating  the  constant  stress  sector  ahead  of  the  crack  from  the  centered  fan  sector.  For 
plane-stress  mode  I,  no  solution  (of  any  type)  seems  yet  to  have  been  found. 

On  the  other  hand,  our  small-a  numerical  solutions  to  the  linear-hardening  problem 
seem  to  follow  a  similar  pattern.  For  the  anti-symmetric  modes,  the  following  features 
appear  as  a  — » 0:  the  strength  of  the  singularity  vanishes  as  s  -  -c  a1/2  (c  >  0),  the  angular 
variations  of  the  stress  and  velocity  fields  seem  to  approach  limits  which  agree  with  the 
corresponding  perfectly-plastic  variations,  but  the  radial  dependence  of  the  velocities  (  es/s  ) 
does  not  approach  the  corresponding  perfectly-plastic  variation  ( lnr  ),  whereas  that  of  the 
stresses  does  in  a  pointwise  sense  ( Is  — ►  1 ).  For  plane-strain  mode  I,  s  does  not  approach 
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zero,  and  neither  do  the  angular  variations  of  the  velocities  appear  to  approach  a  limit  as  a 
approaches  zero,  even  though  the  angular  variations  of  the  stresses  appear  to  approach  a 
limit  that  is  reminiscent  of  the  Prandtl  field.  For  plane-stress  mode  I,  s  does  seem  to 
approach  zero  as  a  0,  and  the  angular  variations  of  the  fields  also  seem  to  approach  a 
limit,  but  we  were  unsuccessful  in  relating  this  apparent  limit  to  a  perfectly-plastic  solution 
with  a  similar  configuration. 

In  this  section  we  will  first  attempt  to  corroborate  the  observations  from  our 
numerical  results  that  as  a  — >  0,  s  — ca1/2,  and  the  angular  variations  of  the  fields  of  our 
small-oc  solutions  approach  the  corresponding  perfectly-plastic  variations.  We  will  second 
try  to  understand  the  breakdown  as  a  — »  0  of  the  radial  dependence  of  our  power 
singularity  solution.  For  simplicity  we  will  initially  address  these  two  questions  in  the 
context  of  the  mode  III  problem,  and  then  we  will  try  to  extend  our  conclusions  to  the 
mode  II  problems.  Finally,  we  will  have  some  comments  about  the  mode  I  problems. 

7.1  Small-a  eigenvalue  and  angular  variations  in  mode  111 

In  Section  2  we  solved  numerically  the  eigenvalue  problem  given  by  the  third  order 
system  (2.5)  subject  to  consistency  conditions  that  determined  the  value  of  a,  and  to  the 
boundary  conditions  (2.9)  and  (2.10).  It  yielded  as  eigenvalue  the  strength  of  the 
singularity,  and  as  eigenfunctions  the  angular  variations  of  the  power-singularity  solution 
(2.4)  to  the  linear-hardening  problem.  The  problem  was  solved  for  several  values  of  a 
ranging  between  zero  and  unity.  The  small-a  results  indicate  that  s  — ca1/2,  and  that  the 
angular  variations  of  the  fields  approach  the  corresponding  perfectly-plastic  variations.  We 
would  like  to  demonstrate  that  this  is  indeed  the  case.  In  order  to  accomplish  this  goal 
rigorously,  we  would  need  to  produce  an  analytical  expression  for  the  small-a  solution, 
and  then  take  its  limit  as  a  — » 0.  However,  we  were  not  able  to  find  such  small-a  solution, 
but  we  did  find  a  solution  to  the  eigenvalue  problem  corresponding  to  the  value  of  a  =  0.  If 
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the  small-a  solution  has  a  unique  limit  as  a  — >  0,  then  it  seems  reasonable  that  this  limit 
should  agree  with  our  a  =  0  solution,  which  we  derive  below. 

Equations  (2.5)  are  in  a  form  that  is  inconvenient  to  work  with  in  the  limit  as  a  -»  0, 


and  they  can  be  rewritten  in  the  following  fashion 

y2’  +  (l+s)yi  =0  (7.1) 

(s/a)  sin9  (y1'-y2)  yi  -  (s2/&)  (cos0  ye2  +  sin0  yiy2)  -  y2y3’  -  syjy3  =  0  (7.2) 
y^'-s  [y2y3+  sine  ^^yi'  +  y^’)]  =  °  (7-3) 

It  is  expected  that  as  a  — » 0,  s  -» 0  so  that  for  a  s  0,  we  take  ssO.  Then  (7.1)  takes  the 
simplified  form 

y2'  +  yi  =°  (7-4) 

In  the  elastic  sector  a  =  1,  and  (7.2)  and  (7.3)  reduce  to 

sin0  (y1'-y2)-y3  =  0  (7.5) 

y3'  =  0  (7.6) 

for  otherwise  we  would  have  to  require  that  y}  =  y2  =  0.  The  solution  of  (7.4)  to  (7.6)  is 
yj  =  A  sin0  +  C  cos0  +  B  (6  sin0  +  cos0  ln|sin0[ ) 
y2  =  A  cos0  -  C  sin0  +  B  (6  cos0  -  sin6  ln|sin0| )  (7.7) 

y3  =  B 


In  the  plastic  sectors  a  =  a.  In  order  to  have  an  appropriate  balance  of  terms  in 
(7.2),  we  need  to  let  either  s/a  or  s2/a  approach  a  constant  non-zero  value  as  a  — >  0. 
Our  numerical  results  support  the  second  possibility:  s2/a  — >  k.  Then  (7.3)  allows  us  to 
identify  two  cases: 

(0  yi  «  0  (7.8a) 

In  this  case  (7.4)  yields 

y2=  1  (7.8b) 

where  we  have  picked  a  normalisation,  and  it  follows  from  (7.2)  and  (7.3)  that 

y3  =  -sin0  (7.8c) 


and  k  =  1. 


Then  (7.2)  yields 

o 

ii 

1 

which  together  with  (7.4)  yields 

yi  =  a  sin0  +  c  cos0 

(7.9a) 

y2  =  a  cos0  -  c  sin0 

(7.9b) 

Also 

y3  =  b 

(7.9c) 

In  order  to  construct  a  solution  to  the  full  problem  we  assume  that  equations  (7.8), 
which  obey  the  boundary  conditions  at  0  =  0,  hold  in  the  interval  0  <  0  <  0j,  equations 

(7.7)  hold  in  0!  <  0  <  02,  and  equations  (7.9)  with  a  =  0  (in  order  to  satisfy  the  boundary 
condition  at  0  =  n)  hold  in  02  <  0  <  7t.  We  then  impose  continuity  of  the  stresses  and 

velocity  at  the  two  boundaries,  the  unloading  condition  which  in  this  case  is  superfluous, 
and  the  reloading  condition  which  yields  |c|  =  1.  The  result,  not  very  surprisingly,  agrees 
with  the  angular  variations  of  the  perfectly-plastic  fields  (see  Appendix  B),  and 
demonstrates  (under  the  presumption  that  the  small-a  solution  has  a  unique  limit)  the 
validity  of  the  original  assertion. 

7.2  Small-a  radial  dependence  in  mode  III 

The  radial  dependence  of  the  velocity  in  our  power-singularity  solution  to  the  linear- 
hardening  problem  suffers  from  the  difficulty  that  it  becomes  unbounded  as  a  (and  hence 
s)  approaches  zero.  We  observe  however  that  for  e_1/lsl  « i  «  1 


Is  =  eslnt  =  1  +  s  Ini  +  ... 

(7.10) 

Is/s  =  1/s  +  lnr  +  ... 

(7.11) 

This  suggests  that  adding  solutions  with  eigenvalues  of  opposite  signs  would  yield  the 
correct  (perfectly-plastic)  radial  dependence  for  both  the  velocity  and  the  stresses  in  the 


limit  as  a  —>  0,  provided  that  the  angular  variations  corresponding  to  both  eigenvalues 
approached  the  same  limit.  Moreover,  it  is  apparent  from  the  discussion  in  the  previous 
section  that  if  a  solution  with  the  specified  balance  of  terms  in  the  limit  (s2/a  -4  constant) 
exists,  the  limiting  angular  variations  are  insensitive  to  the  sign  of  the  eigenvalue. 
However,  because  the  linear-hardening  equations  are  nonlinear,  we  cannot  expect  in 
general  that  the  sum  of  two  such  solutions  with  eigenvalues  of  opposite  sign  will  also  be  a 
solution,  but  in  this  problem  it  happens  that  in  the  limit  as  a  — »  0  the  equations  in  the 
plastic  sectors  become  linear  (in  the  elastic  sector  they  are  always  linear)  so  that  the  sum 
yields  in  the  limit  a  solution  to  the  linear-hardening  problem  which  is  precisely  the 
perfectly-plastic  solution.  In  order  to  give  some  plausibility  to  the  above  ideas  we  will 
study  the  following  model  problem. 

Consider  the  eigenvalue  problem  specified  by  equations  (2.2)  and  (2.3)  with  £x  =  a, 
subject  to  the  homogeneous  boundary  conditions 


Tr(r,0)  =  v3(r,0)  =  0 


V(r>ei)  =  0 


(7.12) 


where  0!  =  tc/2.  Two  solutions  to  this  problem  are 

T*  =  t0  Is*  ee  v3*  =  -V  (t0/G)  sin  9  ?'/  s*  (7. 13) 

with  s*=  ±a1/2  as  eigenvalues.  Note  that 

Te’*  =  -s*  T0  (V/R)  cos0  Is*-1  (7.14) 

It  follows  that  xe'*(r,9j)  =  0,  and  that  for  0  <  0  <  0lf  Te'*(r,0)  >  0  (<  0)  for  s*  =  -a,/2 
(+a1/2)  so  that  of  the  two  possible  eigenvalues  only  the  negative  one  is  physically 
acceptable,  since  we  have  assumed  plastic  loading. 

Further  note  that  the  velocity  in  each  solution  becomes  unbounded  as  a  — »  0. 


However,  surprisingly,  the  addition  of  the  two  solutions  yields  an  exact  solution  to  the 
(nonlinear)  governing  equations,  and  boundary  conditions.  Thus 


X*  =  x0  E(£v<x  +  rVa)/2]  ee  =  t0  cosh(a1/2  Ini)  ee 
v3*  =  -V  (t0/G)  sin0  (i''0  -  rVa)  /  (2a1/2)  =  -V  (t0/G)  sin0  sinh(a1/2  lnr)  /  a1/2 


(7.15) 


is  another  solution  to  our  model  problem.  This  solution  was  first  obtained  through  a 
different  procedure  by  Dunayevsky  and  Achenbach  (1982).  Note  that 

te‘*  =  -oc1/2  t0  (V/R)  cos0  sinh(a1/2 lnr)  >  £  (7.16) 

satisfies  V*(r>0i)  =  0,  and  Te'*(r,0)  >  0  for  0  <  0  <  0j  so  that  (7.15)  is  also  a  physical 

solution.  In  addition,  this  solution  has  the  desirable  property  that  it  is  well-defined  in  the 
limit  as  a  — » 0 

x*  =  lim  {  tg  cosh(a1/2  lnr)  eg  }  =  T0ee 

a_>0  (7.17) 

v3*  =  -lim  {  V  (t0/G)  sin8  sinh(a1/2  lnr)  /  a1/2  }  =  -V  (x0/G)  sin0  lnr 
a— >0 

We  remark  that  these  are  the  fields  (corresponding  to  a  centered  fan  sector)  which 
would  be  the  solution  of  the  corresponding  perfectly-plastic  model  problem;  this 
demonstrates  that  the  governing  equations  of  a  linear-hardening  problem,  at  least  for  some 
configurations,  admit  solutions  that  approach  the  corresponding  perfectly-plasti olutions 
in  the  limit  as  a  — » 0. 

In  our  model  problem  we  find  that  the  power  singularity  solution  to  the  linear- 
hardening  equations  is  not  general  enough  to  yield  the  full  perfectly-plastic  solution  in  the 
limit,  but  does  yield  the  correct  angular  variations  of  the  eigenfunctions.  By  adding  an 
appropriate  nonsingular  (mathematical)  solution  to  the  power-singularity  solution  we  are 
able  to  obtain  the  full  perfectly-plastic  solution  in  the  limit 

We  are  now  in  the  position  to  make  an  estimate  for  the  range  of  validity  of  the 
power  singularity  in  the  model  problem.  We  pointed  out  in  (7.10)  that  for  «  £  «  1, 
Is  -  1.  Therefore,  only  when  0  <  £  «  e_1/^a  are  we  justified  in  neglecting  £^“  compared  to 
£-^“  Hence,  the  range  of  validity  of  the  power  singularity  solution  (7.13)  to  the  model 
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problem  is  given  by  0  <  r  «  e~l/yIa,  and  we  remark  that  this  domain  vanishes  as  a  ->  0, 
as  observed  by  Dunayevsky  and  Achenbach  (1982). 

The  properties  of  the  model  problem  are  very  suggestive,  and  we  would  expect  that 
at  least  some  of  these  features  would  apply  to  the  full  problem.  Indeed,  it  appears  that  the 
full  problem  may  have  another  (mathematical)  solution  with  a  positive  eigenvalue  that 
would  have  the  same  limit  for  the  angular  variations  of  the  fields  and  the  radial  dependence 
of  the  stresses  as  the  power  singularity  solution  (and  the  perfectly-plastic  solution). 
However,  the  radial  dependence  of  its  velocity  would  have  a  singular  term  in  a  that  in  the 
limit  has  the  opposite  sign  of  the  corresponding  singular  term  in  the  radial  dependence  of 
the  velocity  of  the  power  singularity  solution.  Then  the  sum  of  these  two  solutions,  which 
is  not  a  solution  to  the  linear-hardening  equations  for  general  a,  would  in  the  limit  as  a  — » 
0  become  a  solution  to  the  full  linear-hardening  problem,  and  agree  with  the  solution  to  the 
perfectly-plastic  problem.  It  is  conceivable  that  one  might  be  able  to  construct  a  uniformly 
valid  solution  to  the  small-a  problem  (analogous  to  (7.15)  in  the  model  problem)  that 
would  rigorously  demonstrate  the  above  ideas,  but  such  a  solution,  if  it  exists,  has  yet  to  be 
found.  On  the  other  hand,  one  might  expect  that  the  behavior  of  such  small-a  solution 
might  not  be  very  different  in  character  from  (7.15),  and  thus  expect  a  similar  range  of 
validity  for  our  power  singularity  solutions  to  the  full  linear-hardening  problem. 

7.3  The  mode  11  problems 

As  we  have  already  pointed  out,  similar  observations  from  our  numerical  small-a 
solutions  were  made  for  plane-stress  and  plane-strain  mode  II.  It  seems  plausible  that  the 
same  ideas  developed  for  the  mode  III  problem  may  apply  here,  but  the  technical 
difficulties  involved  are  even  greater.  Based  on  these  comments  we  expect  the  range  of 
validity  of  our  power  singularity  solutions  to  be  roughly  the  same  as  that  for  mode  HI. 
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7.4  The  mode  1  problems 

For  plane-strain  mode  I,  we  have  already  observed  that  the  eigenvalue  of  the 
power-singularity  solution  does  not  approach  zero  as  a  — »  0,  and  also  that  the  angular 
variations  of  the  velocities  do  not  appear  to  approach  a  limit.  Alternatively,  the  perfectly- 
plastic  solution  involves  discontinuities  in  the  velocity,  so  that  it  might  be  conjectured  that 
due  to  the  fact  that  our  numerical  formulation  implicitly  assumes  continuity  of  the  velocity, 
our  formulation  of  the  linear-hardening  problem  may  become  overdetermined  as  a  — » 0, 
and  thus  yield  a  spurious  limit  for  s. 

For  plane-stress  mode  I,  the  small-a  solutions  seem  to  indicate  that  a  limit  is 
approached  with  fully  continuous  fields.  But  when  we  searched  for  a  perfectly-plastic 
solution  with  that  configuration  we  failed.  This  may  be  an  indication  that  the  small-a 
results  may  not  actually  be  approaching  a  valid  limit,  in  accordance  with  the  plane  strain 
case. 

Hence  for  the  symmetric  mode  problems  we  expect  our  power  singularity  solutions 
to  be  of  limited  relevance  for  small  a,  and  their  ranges  of  validity  to  be  even  smaller  than 
those  of  the  anti-symmetric  modes. 

8.  Concluding  Remarks 

It  was  found  that  reloading  on  the  crack  flanks  of  the  growing  crack  in  the  linear- 
hardening  material  occurs  for  all  values  of  a  less  than  a  critical  value  a*  where  a*  <  1.  On 

the  other  hand,  it  was  also  found  that  the  angular  variation  (fixed  r)  of  the  effective 
stress  oe  becomes  unbounded  on  the  crack  faces  for  any  amount  of  plasticity  in  a  linear- 

hardening  material  (0  <  a  <  1).  These  two  observations  seem  to  be  in  contradiction.  To  see 

that  they  are  really  not,  we  recall  that  the  definition  of  reloading  involves  a  fixed  particle 
(x2  =  x2).  Hence,  even  though  the  angular  variation  of  the  effective  stress  of  a  particle 
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becomes  unbounded  as  the  particle  travels  far  into  the  wake  of  the  crack  (  8  —>  n ),  the 

radial  variation  (rs)  of  the  same  stress  vanishes,  and  the  product  of  the  two,  namely,  the 
effective  stress  of  the  particle  approaches  a  finite  level,  cew.  We  expect  oew  to  be  a 

continuous  function  of  a,  and  we  know  that  it  vanishes  for  the  elastic  problem  (a  =  1).  On 
the  other  hand  reloading  does  not  occur  until  cew  exceeds  the  effective  stress,  Og1,  that  the 
particle  had  at  unloading.  Hence,  it  follows  that  a*  <  1. 

It  was  found  that  the  size  of  the  plastic  reloading  sector  is  small  (less  than  0.5°)  for 
all  cases  except  for  plane-strain  mode  I  where  the  width  of  the  reloading  sector  exceeded 
40°  for  small  a.  We  also  note  that  the  width  of  the  reloading  sector  for  the  plane-stress 
mode  I  problem  is  significantly  smaller  (less  that  0.001°)  than  for  all  the  other  cases.  Thus 
it  turns  out  that,  with  the  exception  of  the  plane-strain  mode  I  problem,  the  effect  of 
reloading  on  the  global  properties  of  the  solution  (the  strength  of  the  singularity,  s,  and  the 
size  of  the  primary  plastic  loading  sector,  0t)  is  small,  so  that  the  approximation  of  AH  is  a 

good  one.  On  the  other  hand,  in  the  plane-strain  mode  1  problem  we  found  that  the  effect  of 
reloading  on  s  and  0!  is  sizeable;  however,  although  the  effect  is  consistent  in  the  sense 
that  it  makes  s  weaker  for  small  a,  it  does  not  change  the  character  of  the  solution 
significantly.  In  any  case  the  intuitive  expectation  that  reloading  should  become  relatively 
more  important  as  the  perfectly-plastic  limit  is  approached  was  verified  for  all  cases. 

The  mode  II  results  agree  in  general  terms  with  those  of  the  simpler  mode  III  case. 
They  show  a  similar  functional  dependence  of  s,  8j,  and  02  on  a,  and  a  similar  small-a 
behavior.  On  the  other  hand,  the  mode  II  results  introduce  the  new  feature  that  for  large  a 
a  small  reloading  sector  floating  in  the  middle  of  the  elastic  unloading  sector  is  observed, 
which  allows  for  interesting  particle  loading  histories,  but  this  feature  disappears  for 
smaller  a. 

Some  progress  was  made  in  understanding  the  small-a  limit  of  the  power 

singularity  solutions  to  the  linear-hardening  problem.  For  the  anti-symmetric  modes  (II  and 
III),  we  are  inclined  to  think  that  the  small-a  solutions  approach  limits  which  are  related  to 


the  perfectly-plastic  solutions.  In  fact,  the  angular  variations  of  the  small- a  limits  approach 
those  of  the  perfectly-plastic  solutions.  However,  due  to  the  singular  nature  (in  a)  of  the 
radial  dependence  of  the  velocity  in  these  solutions,  their  range  of  validity  (0  <  i «  e*1/>/a) 
vanishes  as  a  — >  0.  On  the  other  hand,  for  the  symmetric  mode  problems,  it  appears  that 
our  small-oc  solutions  are  completely  unrelated  to  the  perfectly-plastic  solutions. 

Even  though  the  validity  of  our  small-a  results  may  be  restricted,  the  results  for 
moderate  and  large  strain-hardening  should  prove  to  be  applicable  to  appropriate  materials 
(such  as  reinforced  ceramics),  where  these  near-tip  solutions  may  be  embedded  in  finite- 
element  calculations  of  the  full  fields  of  a  steadily  propagating  crack. 
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Appendix  A.  The  governing  systems  of  O.D£’s. 

Anti-plane  strain 

fj  =  ( sin0  Tj  y1  [  (s/2)  cos0  +  sin8  T2  y2  +  y3  ] 
f2  =  -(s+  1)  yi 

f3  =  -s2  sin0  yx  -  s2  cos0  y2  -  s  (fit-1  -1)  T3  y2  +  s  (2_I  -I)  sin0  yj  y2  fi 

where 

T^tl+to-J-Dyj2] 

T2  =  [l  +  (fi-1-l)(s+l)y12] 

T3  =  [  s  cos0  +  (s+1)  sin©  y{  y2  ] 
and 

yi  =  (yi  /  ye) 

Plane  strain. 

fj  =  (1— s)  y2  -  2s2  cos0  [  (1+v)  +  (3/2)  (2~l  -1)  ]  y3  -  2s  sin0  Tj 

f2  =  -yj  -  s2  cos0  [  z5  +  (2_1  -1)  x5  ]  -  s  sin9  'i2 

f3  =  -(1+s)  y4  +  y5 

f4=  (l/sin0)  T3  {  ...  } 

f5  =  -(s+2)  y3 

f6  =  (l/sin0)  T3  {  ...  } 

where 

T!  =  s(Uv)  y4- (3/2)  &-1-1)  y3fe 

T2  =  s  y3  +  v  (f4  -  2y3  +  f6)  -  (fi"1  -1)  x5  fc 

T3=  {[1  +  (2"1- 1)  x42][1  +  (2-!-l)  x62]  ~  [v- (2-1-1)  x4  x6]2}-1 

U  =  yc’ 

x4  =  y4-(l/2)(y5  +  y6)  x5  =  y5  -  (1/2)  (y4  +  y6)  x6  =  y6  -  (1/2)  (y4  +  y5) 

z4  =  y4  -  v  (y5  +  y6)  z5  =  y5  -  v  (y4  +  y6)  z6  =  y6  -  v  (y4  +  y5) 


Plane  stress. 


fx  =  (1-s)  y2  -  2s2  cos9  [  (1+v)  +  (3/2)  (fir1  -1)  ]  y3  -  2s  sin9  Ti 
f2  =  -yj  -  S2  cos0  [  z5  +  (ST1  -1)  x5  3  -  s  sine  T2 
f3=-(l+s)  y4  +  y5 

f4  =  (i/sine)  [i  +  (cr1  -l)  V3’1  { yi  + s  cos0  T3 +sine  T4) 
f5  =  -(s+2)  y3 

where 

Tt  =  s  (1+v)  y4- (3/2)  (fit-i -l)y3fe 
T2  =  s  y3  +  v  (f4  -  2y3)  —  (fit  i  — 1)  *5  fe 
T3  =  [z4  +  (fit-1  -1)  x4] 

T4  =  [(2-sv)  y3  -  (fir1  -1)  (3y3f3  +  x5f5  )  x42] 

U  =  Ve’ 

x4  =  y4-(l/2)y5  x5  =  y5  -  (1/2)  y4 

z4  =  y4-vy5  z5  =  y5-vy4 
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where 


Appendix  B.  The  perfectly-plastic  solution  for  anti-plane  strain 

Centered  fan  plastic  sector  (0  <  8  <  0i) 

*  =  *0ee 

v3  =  -V  (T0/G)  sin0  Ini 

Elastic  unloading  sector  (0j  <  0  <  ©2) 

T  =  T0  [  (C  +  B  ln|sin0| )  ej  +  (A  +  B  0  )  e2  ] 
v3  =  V  (t0/G)  B  lm 

Constant  stress  sector  (02  <  0  <  7l) 

*  =  t0*l 

v3=  V  (T0/G)  B  lnr 

A  =  cos0 j  +  0t  sin0j 
B  =  -sin0j 

C  =  sin0j  [  ln|sin0j|  -  l] 

0j  =  19.7112” 

02  -  179.6334” 


Tables 


Table  2.1.- Mode  HI 


a 

s 

s* 

6i 

0i 

1. 

-0.5 

-0.5 

90.00 

90.0 

0.75 

-0.455 

87.79 

0.67 

-0.437 

86.83 

0.5 

-0.394 

-0.394 

84.39 

84.4 

0.3 

-0.325 

-0.325 

79.81 

79.8 

0.2 

-0.277 

-0.277 

76.09 

76.2 

0.1 

-0.207 

-0.207 

69.82 

70.0 

0.05 

-0.153 

63.98 

0.01 

-0.0733 

-0.0737 

52.88 

52.0 

0.005 

-0.0528 

49.21 

0.001 

-0.0244 

-0.0244 

42.88 

43.0 

0.000001 

-0.00081 

33.88 

179.999 

179.947 

179.912 

179.840 

179.731 


(*)  Results  from  Amazigo  and  Hutchinson  neglecting  reloading. 
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Table  3.1.-  Plane  strain  mode  I  v  =  1/3 


a 

s 

s* 

01 

er 

02 

1. 

-0.5 

-0.5 

88.62 

88.6 

0.75 

-0.478 

91.91 

0.67 

-0.469 

93.65 

0.5 

-0.442 

-0.442 

98.40 

98.4 

0.3 

-0.373 

-0.373 

107.46 

107.4 

0.2 

-0.300 

113.88 

179.68 

0.1 

-0.197 

-0.197 

123.35 

124.6 

173.61 

0.05 

-0.142 

-0.136 

130.61 

137.1 

160.95 

0.01 

-0.0797 

-0.0887 

135.19 

156.8 

145.47 

0.005 

-0.0676 

136.07 

141.87 

0.001 

-0.0560 

137.00 

138.36 

(*)  Results  from  Amazigo  and  Hutchinson  neglecting  reloading. 
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Table  4.1.-  Plane  strain  mode  Ii  v  =  1/3 

I 


a 

s 

e, 

e3 

e4 

1. 

-0.5 

42.22 

128.76 

137.66 

0.75 

-0.464 

41.44 

129.45 

136.94 

0.67 

-0.449 

41.10 

129.99 

136.63 

0.5 

-0.413 

40.22 

132.79 

135.85 

0.3 

-0.349 

38.43 

0.2 

-0.301 

36.89 

0.1 

-0.229 

34.12 

0.05 

-0.170 

31.39 

0.01 

-0.0822 

25.90 

0.005 

-0.0594 

24.03 

0.001 

-0.0275 

20.76 

180.00 

179.97 

179.95 

179.91 


Table  5.1.- Plane  stress  mode  I  v  =  l/2 


a 

s 

s* 

ei 

V 

02 

1. 

-0.5 

-0.5 

79.92 

79.9 

0.75 

-0.468 

-0.468 

80.78 

80.8 

0.67 

-0.455 

80.90 

0.5 

-0.420 

80.76 

0.3 

-0.357 

79.34 

0.2 

-0.310 

78.58 

0.1 

-0.237 

-0.237 

73.65 

73.6 

0.05 

-0.178 

-0.178 

69.53 

69.6 

0.01 

-0.0863 

-0.0864 

61.09 

61.1 

180.000 

0.005 

-0.0623 

-0.0624 

58.23 

58.3 

179.999 

0.001 

-0.0287 

53.20 

179.999 

0.0001 

-0.00925 

49.00 

180.000 

(*)  Results  from  Amazigo  and  Hutchinson  neglecting  reloading. 


16.35 

102.58 

134.13 

15.81 

101.61 

132.83 

15.57 

101.33 

132.28 

14.93 

100.96 

130.91 

13.62 

101.37 

128.46 

12.44 

102.62 

126.57 

106.56 

123.50 

58.01 

114.84 

120.34 

180.00 

179.99 

179.91 

179.87 

179.79 
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ft  < o <  1  a  < a < ft 

(d) 

Figure  2.2— <i)  Loading  history  of  particle  P:  it  loads  plastically  (L)  ahead  of  the  crack, 
unloads  elastically  (U)  at  I.  and  reloads  plastically  (L)  at  2.  Loading  history  as  a  function 
of  a  for  modes  III  and  I  (b),  plane-strain  mode  II  (c),  and  plane-stress  mode  □  (d). 
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Figure  2.3 — Strength  of  the  singularity  in  anti-plane  shear. 


Figure  2.4 — Particle  (»2  fixed)  effective  stress  distribution  in  anti-plane  shear  for  large, 
moderate,  and  small  strain-hardening,  normalised  such  that  t,(r,6|)  •  1. 
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Figure  2.Jb— Angultr  i tress  and  velocity  distribution  in  anti-plane  shear  for  moderate 
strain-hardening,  normalised  such  that  I. 
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Figure  2  5c — Angular  unit  and  velocity  distribution  in  anti-plane  shear  for  small  strain- 
hardening,  normalised  such  that  t,(8,)«  1. 
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Figure  2  6 —  Angular  stress  and  velocity  distribution  in  ami-plane  shear  for  perfect- 
plasticity.  normalised  such  that  l„  •  I. 
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Figure  3.1 — Strength  of  the  singularity  in  plane  strain  mode  I  (v  -  1/3). 


Figure  3.2 — Particle  effective  stress  distribution  in  plane  strain  mode  I  for  large,  moderate, 
and  small  strain-hardening,  normalised  such  that  o,(r,6|)  •  I . 
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Figure  3.3a — Angular  stress  distribution  in  plane  strain  mode  I  tor  large  strain-hardening, 
normalised  such  that  0,(6,)  -  1. 


Figure  3. 3h— Angular  stress  distribution  in  plane  strain  mode  1  for  moderate  strain- 
hardening,  normalised  such  that  o,(8,l  ■  I. 
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Figure  3.3c — Angular  sQreas  distribution  in  plane  strain  mode  I  for  small  strain-hardening, 
normalised  such  that  0,(6,)  -  I. 


Figure  3.4 — Angular  stress  distribution  in  plane  strain  mode  I  for  small  strain  hardening, 
normalised  such  that  0,(6,)  •  1. 
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Figure  3.3 — Angular  velocity  distribution  in  plane  attain  mode  I  for  large,  moderate,  and 
small  strain-hardening,  normalised  such  that  ot(8|)  •  I. 
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Figure  4.S — Angular  velocity  distribution  in  plane  strain  mode  U  for  large,  moderate,  and 
small  strain-hardening,  normalised  such  that  o.fOj)  -  I. 


Singularity, 


Figure  5. 1 — Strength  of  the  singularity  in  plane  attest  mode  I  (v  -  1/2). 


Figure  5.2 — Panicle  effective  stress  distribution  in  plane  stress  mode  I  for  targe,  moderate, 
and  small  strain-hardening,  normalised  such  that  o,(r,6|)  •  I. 
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Figure  S.S — Angular  velocity  distribution  in  plane  strew  mode  I  for  large,  moderate,  and 
small  stnin-hardening,  normalised  such  that  Q.(8i)  •  1. 
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Figure  6.1 — Strength  of  the  singularity  in  plane  stress  mode  II  (v  »  1/2). 


Figure  6.2 — Panicle  effective  stress  distribution  in  plane  stress  mode  II  for  large, 
moderate,  and  small  strain-hardening,  normalised  such  that  Ot(r,6|)  -  I. 
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Figure  6  5 — Angular  velocity  distribution  in  plane  stress  mode  II  for  large,  moderate,  and 
small  strain-hardening,  normalised  such  that  Ot(8t)  -  I. 


Figure  6.6— Angular  velocity  distribution  in  plane  stress  mode  II  for  perfectly-piasticity, 
normalised  such  that  o,  «  1. 


